A note on quantum chaology and gamma approximations to 
eigenvalue spacings for infinite random matrices 
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■ Abstract 
X> . 

i Quantum counterparts of certain simple classical systems can exhibit chaotic behaviour through 

Pj^ ■ the statistics of their energy levels and the irregular spectra of chaotic systems are modelled by 

eigenvalues of infinite random matrices. We use known bounds on the distribution function for 
eigenvalue spacings for the Gaussian orthogonal ensemble (GOE) of infinite random real symmet- 
ric matrices and show that gamma distributions, which have an important uniqueness property, 
can yield an approximation to the GOE distribution. That has the advantage that then both 
, chaotic and non chaotic cases fit in the information geometric framework of the manifold of 

gamma distributions, which has been the subject of recent work on neighbourhoods of random- 
^ fin ness for general stochastic systems. Additionally, gamma distributions give approximations, to 

eigenvalue spacings for the Gaussian unitary ensemble (GUE) of infinite random hermitian ma- 
trices and for the Gaussian symplectic ensemble (GSE) of infinite random hermitian matrices 
with real quaternionic elements, except near the origin. Gamma distributions do not precisely 
model the various analytic systems discussed here, but some features may be useful in studies of 
qualitative generic properties in applications to data from real systems which manifestly seem to 
exhibit behaviour reminiscent of near-random processes. 

Keywords: Random matrices, GOE, GUE, GSE, quantum chaotic, eigenvalue spac- 
' ing, statistics, gamma distribution, randomness, information geometry 
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P<j ! 1 Introduction 

O . 

Berry introduced the term quantum chaology in his 1987 Bakerian Lecture [7] as the study of semi- 
classical but non-classical behaviour of systems whose classical motion exhibits chaos. He illustrated it 
with the statistics of energy levels, following his earlier work with Tabor [5] and related developments 
from the study of a range of systems. In the regular spectrum of a bound system with n > 2 degrees 
of freedom and n constants of motion, the energy levels are labelled by n quantum numbers, but the 
quantum numbers of nearby energy levels may be very different. In the case of an irregular spec- 
trum, such as for an ergodic system where only energy is conserved, we cannot use quantum number 
labelling. This prompted the use of energy level spacing distributions to allow comparisons among 
different spectra [E]. It was known, eg from the work of Porter [21], that the spacings between energy 
levels of complex nuclei and atoms with n large are modelled by the spacings of eigenvalues of random 
matrices and that the Wigner distribution [27] gives a very good fit. It turns out that the spacing 
distributions for generic regular systems are negative exponential, that is random; but for irregular 
systems the distributions are skew and unimodal, at the scale of the mean spacing. Mehta [18) pro- 
vides a detailed discussion of the numerical experiments and functional approximations to the energy 
level spacing statistics, Alt et al pQ compare eigenvalues from numerical analysis and from microwave 
resonator experiments, also eg. Bohigas et al [10] and Soshnikov [24] confirm certain universality 
properties. Also, Miller [l!j] provides much detail on a range of related number theoretic properties, 
including random matrix theory links with L-functions. Forrester's online book [15] gives a wealth of 
analytic detail on the mathematics and physics of eigenvalues of infinite random matrices for the three 
ensembles of particular interest: Gaussian orthogonal (GOE), unitary (GUE) and symplectic (GSE), 



o 



- I - 1 



1 



Quantum chaology and gamma manifold approximations 



being the real, complex and quaternionic cases, respectively. The review by Deift p]2j illustrates how 
random matrix theory has significant links to a wide range of mathematical problems in the theory 
of functions as well as to mathematical physics. From 15], we have definitions for the three cases of 
interest: 

GOE: A random real symmetric n x n matrix belongs to the Gaussian orthogonal ensemble (GOE) 
if the diagonal and upper triangular elements are independent random variables with Gaussian 
distributions of zero mean and standard deviation 1 for the diagonal and -j= for the upper 
triangular elements. 

GUE: A random hermitian n x n matrix belongs to the Gaussian unitary ensemble (GUE) if the 
diagonal elements rrijj (which must be real) and the upper triangular elements rrijk — Uj k + 
ivjk are independent random variables with Gaussian distributions of zero mean and standard 
deviation -4= for the rrijj and | for each of the Ujk and Vjk. 

GSE: A random hermitian n x n matrix with real quaternionic elements belongs to the Gaussian 
symplectic ensemble (GSE) if the diagonal elements Zjj (which must be real) are independent 
with Gaussian distribution of zero mean and standard deviation ^ and the upper triangular 
elements Zjk — Uj k + ivj k and Wjk = Uj k + iv'j k are independent random variables with Gaussian 
distributions of zero mean and standard deviation for each of the Ujk, u'j k , v'- k and Vj k . 

Then the matrices in these ensembles are respectively invariant under the appropriate orthogonal, 
unitary and symmetric transformation groups, and moreover in each case the joint density function 
of all independent elements is controlled by the trace of the matrices and is of form [T5] 



where A n is a normalizing factor. Barndorff-Nielsen et al [4] give some background mathematical 
statistics on the more general problem of quantum information and quantum statistical inference, 
including reference to random matrices. 

Here we show that gamma distributions provide approximations to eigenvalue spacing distributions 
for the GOE distribution comparable to the Wigner distribution at the scale of the mean and for 
the GUE and GSE distributions, except near the origin. That may be useful because the gamma 
distribution has a well- understood and tractable information geometry [3J [UJ as well as the following 
important uniqueness property: 

Theorem 1.1 (Hwang and Hu [16]) For independent positive random variables with a common 
probability density function /, having independence of the sample mean and the sample coefficient of 
variation is equivalent to f being the gamma distribution. 

It is noteworthy also that the non-chaotic case has an exponential distribution of spacings between 
energy levels and that the sum of n independent identical exponential random variables follows a 
gamma distribution and moreover the sum of n independent identical gamma random variables follows 
a gamma distribution; furthermore, the product of gamma distributions is well-approximated by a 
gamma distribution. 

From a different standpoint, Berry and Robnik [5] gave a statistical model using a mixture of energy 
level spacing sequences from exponential and Wigner distributions. Monte Carlo methods were used 
by Caer et al. [TT] to investigate such a mixture. Caer et al. established also the best fit of GOE, 
GUE and GSE unit mean distributions, for spacing s > 0, using the generalized gamma density which 
we can put in the form 



p(X) = A n e 



-\TrX~ 2 



(1.1) 



a(f],uj) S 13 e - h{Ms " for ft w >0 



(1.2) 



where a(/3,u>) 




T((2 + /3)M 
J((l +/?)/«;) 



Then the best fits of (|1.2p had the parameter values [H] 



(j. L.J. Dodson 



6 




12 3 4 



Figure 1: The bounds on the normalized cumulative distribution function of eigenvalue spacings for 
the GOE of random matrices h2. 5}) (dashed), the Wigner surmise i2. 3\) (thin solid) and the unit mean 
gamma distribution fit to the true GOE distribution from Mehta JTSjj Appendix A. 15 (thick solid). 
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and were accurate to within ~ 0.1% of the true distributions from Forrester [T5]. Observe that the 
exponential distribution is recovered by the choice g(s; 0, 1) = e~ s . These distributions are shown in 
Figure [3] along with corresponding fits of the gamma distribution. 

2 Eigenvalues of Random Matrices 

The two classes of spectra are illustrated in two dimensions by bouncing geodesies in plane billiard 
tables: eg in the de-symmetrized 'stadium of Bunimovich' with ergodic chaotic behaviour and irregular 
spectrum on the one hand, and on the other hand in the symmetric annular region between concentric 
circles with non-chaotic behaviour, regular spectrum and random energy spacings [81 [TUI [TBI 17]. 

It turns out that the mean spacing between eigenvalues of infinite symmetric real random matrices — 
the so called Gaussian Orthogonal Ensemble (GOE) — is bounded and therefore it is convenient to 
normalize the distribution to have unit mean; also, the same is true for the GUE and GSE cases. 
Barnett [3J provides a numerical tabulation of the first 1,276,900 eigenvalues. In fact, Wigner [35] 
[2l)l [2"7] had already surmised that the cumulative probability distribution function at the scale of the 
mean spacing should be of the form: 

W(s) = l-e~*r (2.3) 
which has unit mean and variance rs 0.273 with probability density function 

77 7TS 2 

w(s) =-se 4 . (2.4) 

Remarkably, Wigner 's surmise gave an extremely good fit with numerical computation of the true 
GOE distribution, cf. Mehta [T5] Appendix A. 15, and with a variety of observed data from atomic 
and nuclear systems [271 151 171 [Tg] . 

From Mehta |18j p 171, we have bounds on the cumulative probability distribution function P for the 



4 



Quantum chaology and gamma manifold approximations 




Figure 2: Probabilty density function for the unit mean gamma distribution fit (thick solid) to the true 
GOE distribution from Mehta \18\/ Appendix A. 15 (points), the Wigner surmised density (thin 
solid) and the probability densities for the bounds \2.1U\) (dashed) for the distribution of normalized 
spacings between eigenvalues for infinite symmetric real random matrices. 



spacings between eigenvalues of infinite symmetric real random matrices: 

L( 8 ) = 1 - e-ra^ 2 < P(s) < U{s) = 1 - e^*'"' (l-^±-\. (2.5) 



Here the lower bound L has mean w 1.13 and variance 4( -^ 7r - ) w 0.348, and the upper bound U 
has mean « 0.940 and variance 96 ~ 7r 2 / 7r « 0.197. 

The family of probability density functions for gamma distributions with dispersion parameter n > 
and mean n/v > for positive random variable s is given by 

f(s;v, K ) = v«^—e- SI/ , for^, K >0 (2.6) 

1 (K) 

with variance Then the subset having unit mean is given by 

/(s;«,k) = K K ^—e- SK , for K >0 (2.7) 

1 (K) 

with variance ^ . These parameters z/, k are called natural parameters because they admit presentation 
of the family (|2.6p as an exponential family 2 and thereby provide an associated natural affine 
immersion in K 3 1141 



.3 . f " 

K 



V 

k I . (2.8) 



iogr(«) — Kiogv 



This affine immersion was used [5] to present tubular neighbourhoods of the 1-dimensional subspace 
consisting of exponential distributions (k = 1), so giving neighbourhoods of random processes. The 
maximum entropy case has k = 1 and corresponds to an underlying Poisson random event process 
and so models spacings in the spectra for non-chaotic systems; for k > 1 the distributions are skew 
unimodular. The unit mean gamma distribution fit to the true GOE distribution from Mehta [18j has 
variance w 0.379 and hence k w 2.42. 
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Figure 3: Probability density functions for the unit mean gamma distributions (dashed) and generalized 
gamma distribution (solid) fits to the true variances for left to right the GOE , GUE and GSE cases. 
The two types coincide in the exponential case, e~ s , shown dotted. 



In fact, k is a geodesic coordinate in the Riemannian 2-manifold of gamma distributions with Fisher 
information metric [3]; arc length along this coordinate from k = a to n = b is given by 



d 2 io g (r( K )) _ i dK 



dn 2 



(2.9) 



Plotted in Figure Q] are the cumulative distributions for the bounds (|2.5p (dashed) , the gamma dis- 
tribution (thick solid) fit to the true GOE distribution with unit mean, and the Wigner surmised 
distribution (|2.3|) (thin solid) . The corresponding probability density functions are in Figure [U unit 
mean gamma distribution fit (thick solid) to the true GOE distribution from Mehta [TH] Appendix 
A. 15 (points), the Wigner surmised density ()2.4j) (thin solid) and the probability densities for the 
bounds (|2.10p (dashed), respectively, 



K s ) = ye 



u(s) = 



tt 2 s(64 



384 



(2.10) 



3 Deviations 

Berry has pointed out [B] that the behaviour near the origin is an important feature of the ensemble 
statistics of these matrices and in particular the GOE distribution is linear near the origin, as is the 
Wigner distribution. Moreover, for the unitary ensemble (GUE) of complex hcrmitian matrices, near 
the origin, the behaviour is ~ s 2 and for the symplectic ensemble (GSM, representing half-integer spin 
particles with time-reversal symmetric interactions) it is ~ s 4 . 

From (|2.7p we see that at unit mean the gamma density behaves like s K_1 near the origin, so linear 
behaviour would require n — 2 which gives a variance of - = | whereas the GOE fitted gamma 
distribution has k ~ 2.42 and hence variance « 0.379. This may be compared with variances for the 
lower bound I, « 0.348, the upper bound u, 96 9 J n » 0.197, and the Wi gner distribution w, 

w 0.273. The gamma distributions fitted to the lower and upper bounding distributions have, 
respectively, kl — ~ 3.660 and k\j — 9 q^ 25 ^ ~ 2.826. Figure [3] shows the probability density 
functions for the unit mean gamma distributions (dashed) and generalized gamma distribution (solid) 
fits to the true variances for left to right the GOE , GUE and GSE cases; the two types coincide 
in the exponential case, e _s , shown dotted. The major differences are in the behaviour near the 
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Figure 4: The unit mean gamma distributions corresponding to the random (non-chaotic) case, k — 
v = 1 and those with exponent k = v = 2.420, 4.247, 9.606 for the best fits to the true variances of 
the spacing distributions for the GOE, GUE and GSE cases, as points on the affine immersion in R 3 
of the 2-manifold of gamma distributions. 



origin. Figure 0] shows unit mean gamma distributions with n = v = 2.420, 4.247, 9.606 for the best 
fits to the true variances of the spacing distributions for the GOE, GUE and GSE cases, as points 
on the affine immersion in R 3 of the 2-manifold of gamma distributions, cf. [Hj. The information 
metric provides information distances on the gamma manifold and so could be used for comparison 
of real data on eigenvalue spacings if fitted to gamma distributions; that may allow identification of 
qualitative properties and represent trajectories during structural changes of systems. 

The author is indebted to Rudnick [55] for pointing out that the GUE eigenvalue spacing distribution 
is rather closely followed by the distribution of zeros for the Riemann zeta function; actually, Hilbert 
had conjectured this, as mentioned along with a variety of other probabilistic aspects of number 
theory by Schroeder [23] . This can be seen in Figure [5] which shows with unit mean the probability 
distribution for spacings among the first 2,001,052 zeros from the tabulation of Odlyzko [50] (large 
points), that for the true GUE distribution from the tabulation of Mehta [18J Appendix A. 15 (medium 
points) and the gamma fit to the true GUE (small points), which has k ~ 4.247. The grand mean 
spacing between zeros from the data was « 0.566, the coefficient of variation ss 0.422 and variance 
w 0.178. 

Table 1 shows the effect of location on the statistical data for spacings in the first ten consecutive 
blocks of 200,000 zeros of the Riemann zeta function normalized with unit grand mean; Table 2 
shows the effect of sample size. For gamma distributions we expect the coefficient of variation to be 
independent of sample size and location, by Theorem ll.il 

Remark The gamma distribution provides approximations to the true distributions for the spacings 
between eigenvalues of infinite random matrices for the GOE, GUE and the GSE cases. However, it 
is clear that gamma distributions do not precisely model the analytic systems discussed here, and do 
not give correct asymptotic behaviour at the origin, as is evident from the results of Caer et al. [11] 
who obtained excellent approximations for GOE, GUE and GSE distributions using the generalized 
gamma distribution (|1.2p . The differences may be seen in Figure [3] which shows the unit mean 
distributions for gamma (dashed) and generalized gamma [llj (solid) fits to the true variances for the 
Poisson, GOE, GUE and GSE ensembles. 
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Figure 5: Probability plot with unit mean for the spacings between the first 2,001,052 zeros of the 
Riemann zeta function from the tabulation of Odlyzko \20f/ (large points), that for the true GUE 
distribution from the tabulation of Mehta [18] Appendix A. 15 (medium points) and the gamma fit to 
the true GUE (small points). 
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Table 1: Effect of location: Statistical data for spacings in the first ten consecutive blocks of 200,000 ze- 
ros of the Riemann zeta function normalized with unit grand mean from the tabulation of Odlyzko \20}/ . 
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Table 2: Effect of sample size: Statistical data for spacings in ten blocks of increasing size 
200, 000m, to = 1,2,..., 10, for the first 2,000,000 zeros of the Riemann zeta function, normalized 
with unit grand mean, from the tabulation of Odlyzko \20f . 

The generalized gamma distributions do not have a tractable information geometry and so some 
features of the gamma distribution approximations may be useful in studies of qualitative generic 
properties in applications to data from real systems. It would be interesting to investigate the extent 
to which data from real atomic and nuclear systems has generally the qualitative property that the 
sample coefficient of variation is independent of the mean. That, by Theorem 11.1) is an information- 
theoretic distinguishing property of the gamma distribution. 

It would be interesting also to know if there is a number-theoretic property that corresponds to 
the apparently similar qualitative behaviour of the spacings of zeros of the Riemann zeta function, 
Tables 12.1, 12.2. Since the non-chaotic case has an exponential distribution of spacings between 
energy levels and the sum of n independent identical exponential random variables follows a gamma 
distribution and moreover the sum of n independent identical gamma random variables follows a 
gamma distribution, a further analytic development would be to calculate the eigenvalue distributions 
for gamma or loggamma-distributed matrix ensembles. Information geometrically, the Riemannian 
manifolds of gamma and loggamma families are isometric, but the loggamma random variables have 
bounded domain and their distributions contain the uniform distribution, which may be important in 
modelling some real physical processes. 
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